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ABSTRACT 

We study the disk emission component hidden in the single-peaked Broad Emission 
Lines (BELs) of Active Galactic Nuclei (AGN). We compare the observed broad lines 
from a sample of 90 Seyfert 1 spectra taken from the Sloan Digital Sky Survey with 
simulated line profiles. We consider a two-component Broad Line Region (BLR) model 
where an accretion disk and a surrounding non-disk region with isotropic cloud veloc- 
ities generate the simulated BEL profiles. The analysis is mainly based in measure- 
ments of the full widths (at 10%, 20% and 30% of the maximum intensity) and of the 
asymmetries of the line profiles. Comparing these parameters for the simulated and 
observed Ha broad lines, we found that the hidden disk emission may be present in 
BELs even if the characteristic of two peaked line profiles is absent. For the available 
sample of objects (Seyfert 1 galaxies with single-peaked BELs), our study indicates 
that, in the case of the hidden disk emission in single peaked broad line profiles, 
the disk inclination tends to be small (mostly i < 25°) and that the contribution 
of the disk emission to the total flux should be smaller than the contribution of the 
surrounding region. 
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1 INTRODUCTION 

Modeling of the double-peaked Balmer lines has been 

kinematics of the Broad 



Line Region (BLR) (see e.g. 


Perez et al. 19881; IChen et all 


198C 


Chen & HalDern 1989 


; Eraclcous & Haloern 1994 


2002 


Rodriguez- Ardila et al.l 


1996; Storchi-BcrEmann et al.l 


1997 


Livio & Xulll997'; 'Ho ct al. 2000; Shields et al. 2000l; 


Strateva et al. |2003: Storchi-Bcremann et al. 2003albl). 



However, only a small fraction of Active Galactic Nuclei - 
AGN (3%-5%) shows clearly d ouble peaked B road Emission 
Lines (BELs) in their spectra (|Strateva et al.| [2003l. 

According t o th e standard unification model 
l|Urrv fc Padovanil Il995l ) one can expect an accretion 
disk around a supermassive black hole in all AGN. The 
majority of AGN with BELs, have only single peaked 
lines, but this does not necessarily indicate that the 
contribution of the disk emission to the BELs profiles is 
negligible. It is well known that a face-on disk also emits 
single peaked broad lines ( see e. g. Chen fc HalpernI Il989l; 
iDumont fc CoUin-SouffrinI I199G| ; Kollatschnv fc Bischofj 



|2002| ; iKoUatschnvl |2003| ). Moreover, a Keplerian disk 
of arbitrary inclination with presence of a disk wind 
can also produce sin gle-peaked broad emission lines 
iMurrav fc Chianel[l995l ). 

In spite that most of the BELs are single peaked, 
there are other evidences like the detection of asymme- 
tries and substructure (shoulders or bumps, for instance) 
in the line profiles that indicate the presence of a disk (or 
disk-like) emission JPopovic et al.l [20021 ; iKollatschnvl l2003l : 
IShapovalova et al] 12004 ). Also, the study of the accretion 
rates in AGN supports the presence o f a standard opt i- 
cally thick and geometrically thin disk l|Wang et al.ll2003l ). 
Moreover, the spectropolarimetric observations gave an ev- 
idence for the disk -like emission (rotational motion, see e.g. 
ISmith et al1l2005l ) 

To explain the complex morphology of the observed 
BELs shapes, different geo metrical models have been dis- 
cussed (see in more details ISulentic et alll2000l ). In some 
cases the BELs profiles can be explained only if two or more 
kinematically different emission regions are considered (see 
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e.g. 


Romano et al. 19961: PoDovic et alj 2001 , 20021. 


2003 
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. 2008. 


200g|; 


Bon et al. 


20061: Bonll2008l:lllic et all 


2006 


Collin et al 


|200€ 


: IHu et al. 


I2OO8I'). In Darticular. the 


exis- 
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tence of a Very Broad Line Region ( VLBR) with random ve- 
locities at 5000-6000 km/s within an Intermediate Line Re- 
gion (ILR) has also been considered to explain the observed 
BELs profiles (jCorbin fc Boroson|[l996l : ISulentic et allbOOOl : 
IHu et al. ''2008Y 

In this paper we study the presence of the hidden disk 
emission in objects which show only one dominant peak in 
their broad emission line profiles. To do that we consider 
that the BLR has two kinematic components; an accretion 
disk and a surrounding non-disk region. With this model 
we compute emission line profiles for different values of the 
model parameters. Then we compare the simulated profiles 
(specifically, their widths and asymmetries) with the obser- 
vational data. 

The aim of this paper is to discuss possibility that the 
disk geometry, at least partly, affects the complex BEL pro- 
files, i.e. to try to constrain the BLR geometry that is im- 
portant for estimates of the AGN black ho le masses and 
accretion rates (see e.g. |La Mura et al]|20Q9l ). 

The paper is organized as follow: in §2 we describe the 
two component model of the BLR and perform the numerical 
simulation. In §3 we compare the simulations with available 
data. In §4 we discuss our results and in §5 we outline our 
conclusions. 




Figure 1. The scheme of the BLR composed from two geomet- 
rically different components: accretion disk and clouds randomly 
distributed around the disk. 



The local broadening (a) and shift izi) of each disk 
element have been taken into account as in lChen fc HalpernI 
(1989), i.e. the 5 function has been replaced by a Gaussian 
function: 



2 NUMERICAL SIMULATIONS BASED ON A 
TWO-COMPONENT MODEL 

2.1 BLR geometry 

In the last years, arguments supporting the presence of disk 
winds show ability to explain a number of observed AGN 
phenomena such as the X-ray and UV absorption, line emis- 
sion, reverberation results, some differences among Seyfert 
and other active objects (like quasars or broad-line radio 
galaxies), and the presence or absence of do ublc-pcaked 
emission-line profiles (se e e.e lMurrav fc Chiang. 1995 . 199^ : 
iProga fc Kallmanll2004h . These results support a model in 
that the BLR is composed from two kinemati cally dis- 
tinct components, a disk and a wind. Recently, IHu et al] 
l|2008l l confirmed that the BLR is probably composed 
from two emission regions, i.e. a VBLR and a ILR com- 
ponent, as it was earlier assumed in several papers (see e.g. 
Corbin fc Boroson|[l99^ : ISulentic et al]|2000l : iPopovic et al] 
2004 etc.) . 

Consequently, we assumed that the BELs can be kine- 
matically divided into two components, one from the VBLR 
(contributing to the wings) and other from the ILR (con- 
tributing to the core). Of course, one can assume different 
geom etries for both emitting regions (see e.g. [Popovic ct al. 
|2004| . and also Appendix B in this paper), but here we will 
assume that the VBLR is coming from an accretion disk, 
and the ILR from an additional region, which surrounds 
the disk, and has an isotropically distributed random veloc- 
ity. Note here that the kinematics of a wind would imply 
radial velocities (logarithmic profile), more than isotropic 
(Gaussian profile) , but we assume Gaussian profile as a first 
approximation. The scheme of the assumed model is pre- 
sented in Fig. 1. 



We express the disk dimension in gravitational radii 
(Rg — GMjf? , G being the gravitational constant, M the 
mass of the central black hole, and c the velocity of light). 

On the other hand, we assume that the additional emis- 
sion region can be described by a surrounding region with 
an isotropic velocity distribution, i.e. the emission line pro- 
file generated by this region can be described by a Gaussian 
function with broadening Ws and shift Zg. Thus, the whole 
line profile can be described by the relation: 

J(A) =/,(A) + /,(A), (2) 

where /d(A), and is (A) are the emissions of the relativistic 
accretion disk and the non-disk region, respectively. 

2.2 Parameters for the disk and surrounding 
region 

As it was earlier noted in iPopovic et al.l (|2004l '). this two 
component model can fit the line profiles of BELs, but is too 
open to constrain the physical parameters. First of all, the 
disk model includes many parameters (the size of the emit- 
ting region, the emissivity and inclination of the disk, the 
velocity dispersion of the emitters in the disk, etc.). There- 
fore, in order to do numerical tests, one needs to introduce 
some constraints and approximations. 

It seems that the parameter a of the Doppler broaden- 
ing of the non-disk region and parameter a corresponding to 
the broadening of the random moti o n in the disk mo del are 
connected (see IPopovic et al.|[20o3 : iBon et al] 120061 ). As a 
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first approximation we assume that the random velocities in 
the disk and in the non-disk region are the same. So here we 
consider a parameter a — 1000 km/s for both, the Doppler 
broadening of the non-disk regions a s well as the a in the 
model of the disk profile (also, see lEracleous &: HalpernI 

On the other hand, we considered a wide range of disk 
parameters but with several constraints: 

i) The disk inclination affects the emission obtained 
from the disk. The observed fiux from the disk (Fd) is pro- 
portional to the disk surface (Sd), as 



Sd ■ cos(z) 



(3) 



where i is the inclination, and S^f / is the effective disk emit- 
ting surface, therefore, one cannot expect a high contribu- 
tion of the disk emission to the total line profile for a near 
edge-on projected disk. 

ii) As far as the plasma inside 100 gravitational radii 
from the central black hole is very hot, one cannot expect 
emission of the low ionized lines in this part of the disk. Then 
we limit the inner r adius to Ri„n > 1 Rp. Consequently, 
the model given by IChen fc Halperi] (|l989l ) can be prop- 
erly used, i.e. it is not necessary to include a full r elativistic 
calculation (as e. g. in ljovanovic fc PopovidboOSh . 

iii) The emissivity of the disk as a function of radius, 
r, is given by e — eor~^. Since the illumination is due to 
a point source radiating isotropically, located at the cen- 
ter of the disk, the fiux in the outer disk at di fferent radii 
should vary as r~ (lEracleous fc HalpernI Hooi l. We note 
here that this is indeed the way how the incident flux varies, 
but not necessarily the way in which lines respond to it 
(see e.g. iDumont fc Collin-Souffrinlll990l : iRokaki fc Boissoiil 
1 19991 : lEracleous fc Halpernll2003l ). However, the power in- 



dex p ~ 3 can be adopted as a reasonable approxima- 
tion at least for H q (Eracleous fc Halgcrn ,2003). Also, 
ICao fc Wan3 (|2006l ') indicate that p is in the range from 2 
to 3. Moreover we simulated the influence of the emissivity 
(p) to the disk line profile (Fig. 2) and found that it only 
slightly affected the normalized disk line profiles. Therefore 
the assumption of p = 3 can be accepted for purposes of this 
work. 

iv) Previous estimatio ns of the d oublc-pcak ed AGN 



19891: Chen & Halpern 


1989; 


Eracleous & Halpern 1994; 


Rodriguez- Ardila et al.iil996: Storchi-Bergmann et al, 19971: 


Ho et al. 2000l: Shields et all 


2000l: 


Strateva et al. 20031: 


Storchi-Bergmann et al.l 


2003a Ibl: 


Eracleous & Halpern 


20031) show that the tvpical dimensions of an accretion 



disk that emits low ionization lines are of the order of 
several thousands Rg. For that reason we did not consider 
dimensions of the disk larger than 100000 Rg. This was an 
important approximation to limit the computing time. 

v) We consider a systemic velocity shift of the non-disk 
region (not greater than ± 3000 km/s) to test the possibility 
of the outflow/inflow. 



In Appendix A we give simulations of different a for non-disk 
region 




).05 -0.04 -0.03 -0.02 -0.01 0.01 0.02 0.03 0.04 0.05 0.06 
V/c 

Figure 2. Simulated disk profiles for different values of the emis- 
sivity: p=2 (solid line) p=3 (dashed line) and p=4 (doted-dashed 
line), for two inclinations i =10 (narrower lines) and 30 degrees 
(broader line profiles). The inner and outer radius are taken as: 
Rinn =500 Rg and i?out=1500 Rg. The maximum intensity is 
scaled to one. 



2.3 Flux ratio, normalized widths and asymmetry 
parameters 

We consider the following parameters to study the simulated 
and observed BELs profiles: 

i) The flux ratio between the disk (Fd) and the non-disk 
region (Fs): 



where 



El 

Fd 



Ftot = Fs+Fd = il + Q)Fd 



(4), 



Using this parameter, the total line proflle (normalized 
to the disk fiuxjf] can be written as: 



Itot{\) _ IdW , „/s(A) 



+ I 



Fs 



(5) 



Fd Fd 

where /(A) is the wavelength dependent intensity. The com 
posite profile is normalized according to, 

/tot (A) 



3(A) 



J max 
^tot 



(6) 



profile. 

ii) For the composite line profile 9(A) we measured full 
widths at 10%, 20%, 30% and 50% of the maximum inten- 
sity, i.e. 11)10%, '^20%, "^30% a-iid w^q%. Then we define co- 
efficients ki {i = 10, 20, 30) normalized to the Full Width 
at Half Maximum (FWHM), as fcio = wio%/'W5o%, feo ~ 
W2o%/w^Q% and fcao = W3o%/w5o%. It is obvious that the 

^ This is taken from technical reasons to simulate different con- 
tributions of the disk and the non-disk component. First we 
normalized both line profiles to their fluxes, and after that we 
rescaled the non-disk component multiplying with Q, then whole 
proflle is given in units of the disk fiux. 
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coefficients ki are functions of the radius 7? and other pa- 
rameters of the disk. Using these normalized widths we can 
compare results from AGN with different random velocities. 

iii) We also measured the asymmetry {Ai) at i =10%, 
20%, 30% of maximum intensity of the modeled and ob- 
served lines as 



FWHM 



(7) 



where W/* and VKf are red and blue half widths at i = 10%, 
20% and 30% of the maximum intensity, respectively. 



2.4 Simulated line profiles 

First of all, we simulated only the disk profiles, taking into 
account different values of the disk parameters. An exten- 
sive discussion ab out possible disk line profiles is given in 
iDumont fc Collin-So uffrin (1990,). In the first instance, the 
relative importance of the disk contribution to the core or to 
the wings depends on the disk inclination. In Fig. 3 we pre- 
sented simulated profiles, corresponding to Rinn = 400 Rg, 
Rout ~ 1200 Rg and 12000 Rg, for different inclinations: 
i = 1°, 10°, 20°, 40° and 60°. As it can be seen in Fig. 3, 
the contribution of the disk to the center of the line or to 
the wings is not so much sensitive to the outer radius, but 
significantly depends on the disk inclination. A face-on disk 
contributes more to the core of the line, while a moderately 
inclined disk (40° > i> 20°) contributes significantly to the 
line wings. For i > 40° the disk emission will strongly affect 
the far wings of the composite profile. 

Another very important parameter is the fiux ratio 
between components, Q. As examples, in Fig. 4 we pre- 
sented five simulations of composite line profiles with values 
Q = 0.3, 0.5, 0.7, 1, 1.5 and 2, where Rir,„ = 400 Rg and 
Rout ~ 1200 Rg, and for different inclinations (z =1, 10, 20, 
40, 60 degrees). 

We found that the presence of the disk emission is diffi- 
cult to detect in the line profile when the contribution of the 
disk is smaller than 30% of the total line emission (Q > 2): 
in the case of a low inclination both the disk and non-disk 
region contributes to the line core and it is very hard to sep- 
arate the disk and nod-disk region. In the case of a highly 
inclined disk, the disk emission spreads in the far wings, and 
could not be resolved from the continuum, especially if the 
observed spectrum is noisy. For the case of dominant disk 
emission {Q < 0.3), if the inclination is low, the line will be 
shifted to the red, and if the inclination is high, two peaks 
or at least shoulders should appear in the composite line 
profile. Consequently, further in the paper we will consider 
only cases where 0.3< Q <2. 

Note here that, although we have considered a relatively 
low random velocity of ^1000 km/s, the lines where the disk 
is dominant (the disk contributes at least 50%) can be very 
broad. The obtained widths are in agreement with the mea- 
sured wid ths of double-peaked lines, that ra nge from several 
thousand (|Eracleous fc HalpernI 1994120031 ) to nearly 40,000 



kms ^ (see e.g. lWang et al. 
Ha width of « 40000 km s 



20051 . SDSS J0942-h0900 has the 



2.5 Results from the line profiles simulations 

From the simulations mentioned above we infer the follow- 
ing results: 

(i) To detect disk emission in a BELs, the fraction of the 
flux emitted by the disk in the total line profile should be 
higher than 30%. A dominant disk {Q < 0.3) will be clearly 
present in the total line profile (peaks or shoulders in the 
line profiles). 

(ii) In the case of a nearly face-on disk {i < 5°), the 
disk emission may contribute to a slight asymmetry towards 
the red (due to the gravitational redshift), but it is hard to 
detect this asymmetry. In the case of an edge-on disk, the 
emission from the disk will contribute to the far wings and 
then it may be difficult to separate it from the continuum. 

(iii) These two parameters, the fiux ratio between com- 
ponents and disk inclination, are crucial for the line shapes 
in the two-component model. 

(iv) In the simulated line profiles the asymmetry was 
mostly Ai > 0. For low inclination (i < 10°), the asymmetry 
weakly depends on Q. 



3 COMPARISON BETWEEN SIMULATED 
AND OBSERVED BEL PROFILES 

3.1 Data sample and measurements 

Th e set of spectra for ou r data sample has been collected 
by iLa Mura et al.l (|2007h from the spectral database of 
the third data release from the Sloan Digital Sky Survey 
(SDSSfl According to the purposes of the work (see also 
iLa Mura et al.ll2007l )). the SDSS database was searched for 
sources corresponding to the following requirements: (i) ob- 
jects were with redshifts z < 0.4, so Ha would be covered 
by the available spectral range; (ii) the Balmer series were 
clearly recognized, at least up to H5 in order to see that in all 
Balmer lines a broad component is present; (iii) and profiles 
were not affected by distortions (bad pixels on the sensors, 
the presence of strong foreground or background sources). 

The preview spectra provided by the database retrieval 
software were manually inspected, looking for objects in bet- 
ter agreement with our requirements, until 115 sources were 
chosen from approximately 600 candidates examined in var- 
ious survey areas. Subsequent inspection of the spectra col- 
lected within the database led to the rejection of 25 objects, 
which were affected by problems that could not be detected 
in the preview analysis. Therefore, our resulting sample in- 
cludes the spectra of 90 various broad- line-emitting AGN, 
corresponding to 15% of the candidates that we examined 
and located in the range. 

The spectra were already corrected for instrumental and 
environmental effects, including sky-emission subtraction 
and correction for telluric absorption and calibration of data 
in physical units of the flux and wavelengt h. Spectra were 
corrected for the Galactic extinction (see iLa Mura et al.l 
I2007D . Also, the cosmological redshift calibration were per- 
formed. Since the interest was to investigate the broad line 
shapes, the subtraction of the narrow components of Ha as 
well as the satellite [Nil] lines were performed. The spectral 
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Figure 3. Simulated disk line profiles for five different disk inclinations (i=l, 10, 20, 40 and 60 degrees, from the higfiest to the lowest 
intensity lines, respectively) for the disk with the fixed inner radius Ri„„ =400 Rg and with outer _Rout=1200 Rg (left) and Rout =12000 
Rg. (right). 




Figure 4. Simulated line profiles emitted by the two-component model for five different inclinations (i=l, 10, 20, 40 and 60 degrees, 
from the narrowest to the broadest line, respectively) for different contributions of the disk to the composite line profiles (as it is written 
in figures). The inner radius of the disk is taken to be 400 Rg, while the outer radius is 1200 Rg. 



reduction (including subtraction of stellar component) and 
the way to obtai n the broad line prof iles are in more de- 
tails explained in lLa Mura et all (l2007l ). We have also used 
FWHM and FWZI measurements from Table 2. of the men- 
tioned paper. 

Previously cleaned broad Ha profile was normalized, 
converted from wavelength to velocity scale and smoothed, 
with Gaussian smoothing, using DIPSO software packag^fl- 
Additionally, we have measured the half (red and blue) 
widths, Wi^ and VF/', at 10%, 20% 30% and 50% of max- 
imum intensity (see Fig. 5), and after that we calculated 
normalized full widths, ki, and asymmetries, Ai (Eq. (7)). 



http:/ /www. starlink.rl.ac.uk 



3.2 Observed vs. simulated line profiles 
parameters 

In Fig. 6 we presented the normalized widths k2o vs. fcio and 
fcao vs. fcio measured for Hq lines of the sample. We also 
plot the values corresponding to simulations with different 
disk inclinations, fiux ratios Q, and fixed disk dimensions of 
Rinn ~ 600 Rg and Rout = 4000 Rg. As one can see in Fig. 6, 
most of the measured points are located within 0.7 < Q < 2, 
and i < 25°. These results do not change significantly if we 
change the inner and outer radii of the disk, where the inner 
radius is not taken to be less then 100 Rg, while the outer 
was not taken to be closer then several hundreds Rg. 

The measured asymmetries of the Ha profiles are pre- 
sented in Fig. 7. Lines in this Figure correspond to the sim- 
ulated models as in Fig. 6. As it can be seen in Fig. 7, some 
objects (around 10 objects from the sample of measured Ha 
lines) have a high negative asymmetry which cannot be ex- 
plained by assumed models without a systemic motion of the 
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1.5 2 2.5 3 3.5 2 2.5 3 3.5 

Figure 6. The measured width ratios (crosses) and simulated values (dashed lines) from the two-component model for the different 
contribution of the disk emission to the total line fiux (Q=0.3, 0.5, 0.7, 1, 1.5 and 2). The inner disk radius is taken to be 600 Rg, outer 
4000 Rg, and different inclinations are considered (solid isolines presented i = 10, 15, 20, 25 and 30 degrees, respectively). 



Wb(1 0%) / 


\ Wr(10%) - 







-8000 -6000 -4000 -2000 2000 4000 6000 8000 
velocity (km/s) 

Figure 5. The broad H« line of SDSSJ1025-f 5140. The horizon- 
tal lines presented the measured widths at 10%, 20%, 30% and 
50% of the maximum intensity. The blue and red half widths have 
been measured. 



non-disk region. To explore the outflow or inflow presence 
(to explain the blue asymmetry of the core component), we 
consider shifts of the Gaussian component in the model. As 
it can be seen in Fig. 7 (bottom) models where a blueshift 
of the Gaussian component of 800 km was taken into 
account (other parameters are not changed), are able to ex- 
plain the measured asymmetries. 

As it can be seen in Fig. 6, the coefficients k are sen- 
sitive to the disk inclination and Q, therefore we will use 
the relationships between these coefficients to give some es- 
timates of the inclination and Q for the sample. Of course, 
we should also take into account the influence of the disk 
dimensions. To do that we compute a grid of models with 
ki — f(Rinn,i, Q) confirming that the changes in the param- 



eters mainly depend on i and Q. We found that for inner 
radius Kinn <200 Kg, the simulated ki do not fit well the 
observations (most of the measured points from the sample 
are out of the grid of models, as it is presented in Fig. 6). 
The same inconsistencies are found for Rout <1000 Hg. 

In order to derive some results about the disk parame- 
ters we assume the following constraints and procedures: 

i) According to the res ults obtained from th e fittin g of 
the double-peaked lines in lEracleous fc HalpernI l|2003h we 
fixed the inner radius at Ri„„=600 Kg, and the outer at 
Ront=4 000 Rg, as the averaged va lues obtained from their 
fittings l|Eracleous fc Halpernllioosl ). 

ii) For each AGN in the sample we estimate i and Q 
using the normalized widths, ki. Specifically, we obtain two 
estimates for i and Q as the values associated to the mea- 
surements of both fcio vs. k2o and fcio vs. feao. In Table 1 we 
present the averaged values and differences between those 
estimates, Ai, and AQ for the case without blue shift of the 
non-disk component.. 

iii) We excluded from the analysis the objects where 
difference between estimated i and Q was huge, in total 5 
objects where the two-component model cannot be applied. 
Also there are 9 objects Ai > 1 degree (presented as full 
triangles in Fig. 8 and in Figs, after). In Fig. 8, we present 
the ki parameters as function of the inclination. As one can 
see, most of the points are well concentrated as a linear 
function of i vs. ki. 

In Fig. 9, we present histograms of the number of AGN 
vs. i and Q. As one can see in Figure, there is a peak at 
i ~ 17° while estimated values are within 10° < i < 25° 
in both cases: without (solid line) and with systematic blue 
shift of the non-disk (see Fig. 7). Also, there is a peak at 
Q ~ 1.25 and most of the points are within 1 < Q < 2, 
showing that the disk emission is typically smaller than the 
emission of the non-disk region. 

One can expect randomly oriented accretion disk in 
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Figure 7. The asymmetry of the Ha line of the sample (crosses) and simulated one from the two-component model. Different contributions 
of the disk emission to the total line flux (Q) are taken into account, and plotted with a dashed lines (Q=0.3, 0.5, 0.7, 1, 1.5 and 2). The 
inner disk radius is taken to be 600 Rg and outer 4000 Rg. Different inclinations are considered, taking values from (full squares), 10 
(open circles), 20 (full circles) and 30 (open triangles) degrees. On two panels up we presented the models without systematic shift of 
the non disk component, while on two panels bottom, a systematic blue shift of 800 km/s was taken into account. 




10 15 20 

inclination (degree) 




Figure 9. Histograms of the inclination (left) and Q = Fs/F^^gk (right) for the sample in the case with (dashed line) and without (solid 
line) systematic blue shift of the non-disk component (see Fig. 7) . The cases where discrepancy in the inclination were more than one 
degree, were not taJsen into account. 
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Halpern 1994), where is assumed that in these objects Q < 0.1. 
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Figure 8. Inclination vs. feio (the first), k20 (the second) and 
feso (the third). The points where description were more than 
one degree (mainly above the other points) are denoted with full 
triangles. 



AGN, but we obtained low inclined disk. Such small inclina- 
tion range (10° < i < 25°), may be expected since a highly 
inclined disk has a smaller brightness than surrounding non- 
disk region. Therefore one can expect a weak disk emission 
in far wings which cannot be detected. 



4 DISCUSSION 

In the previous sections we have made a grid of two- 
component models, aiming to search for the hidden disk 
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Figure 11. Inclination (top) and Q (bottom) as function of 
the FWZI/FWHM. With full points the measurements of the 
disk parameters and widths are taJien from Eraclous & Halpern 
(1994,2003), and open circles are from the sample (measurements 
of the FWHM and FWZI arc taken from La Mura et al. 2007). 
The vertical dashed line shows a rough border between disk emis- 
sion lines and ones where the disk may be superposed with the 
emission of an additional region. 



emission. We found that the more significant parameters 
in the generation of the emission line profiles are the incli- 
nation and the flux ratio between components, Q. On the 
other hand, wo have also explored the influence of the inner 
radius. Comparing a grid of simulated line profiles with 
the data of the 90 AGN sample from SDSS, we found that, 
if the disk emission is present, the inner radius should not 
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be smaller than 200 Kg (it seems to be in the range 1000 
Rg > Ri„„ > 200 %fl 

Fixing the inner and outer radius to an averaged value, 
obtained from the st udy of BELs with two-peaked lines 
IIEracleous fc HalpernI [20031. we estimated the values of i 
and Q for the 90 AGN sample from SDSS. According to Ta- 
ble 1 (where we give data for the case without systematic 
blue shift of the non-disk component) and the histogram in 
Figure 9, the two component model associates significant 
disk emission to practically all objects. 

In Figure 10 we present Q vs i, with well known double 
peaked AGN, assuming that the contribution of the non-disk 
region is smaller than 10%. Note here, that even in double- 
peaked emitters, there is often a residual, so called 'classical 
broad hne' c omponent left over after th e disk fit (as e.g. in 
3C390.3, see lEracleous fc Halpern|[2003h . In this case, there 
is an indication that the linear regression may be present 
in Q vs. i. This may be caused by the disk brightness, i.e. 
with higher inclinations the disk emission decreases (Eq. 2) 
and in this case, the disk emission can be detected if the 
non-disk emission is negligible. 

As an additional tes t, we have used meas urements of 
the FWHM and FWZI (|La Mura et all |2007| ) for the 90 
AGN sample from SPSS and fo r double peaked lines by 
lEracleous fc HalpernI l|l994l2003h to plot the inclination, i, 
vs. FWZI/FWHM. As it can be seen in Fig. 11, there 
is a clear separation between single and double-peaked lines 
(except for the points for which we estimate Q and i with 
high uncertainty, presented as full triangles). This also indi- 
cates that a high inclined disk emission can be detected if it 
is more dominant than the non-disk component. Note here 
that the lack of AGN population with Q < 1 (see Fig. 11) 
is probably caused by selection effects, since we selected a 
sample where single-peaked profile is dominant. 

The problem of a low- inclined disk (i < 25°), that we 
obtained from this sample of AGN, still remains. The restric- 
tion on inclination is problematic for disk models, since we 
may expect random orientations - at least within the range 
allowed by the torus opening angle (given that the accretion 
disk and torus are co-axial, see Fig. 12) . Then high-inclined 
disks may be obscured by the torus, but in this case one can 
expect the cut-off in the inclination around 45°-50°. One 
explanation may be that we have situation as schematically 
presented in Fig. 12, where the disk is co-axial with the 
tours, and there are three cases depending on the line-of- 
sight of an observer; 1) line-of-sight is throughout the torus, 
where the disk and non-disk regions are obscured and one 
can detect only narrow lines, 2) as inclination stay smaller, 
the non-disk contribution stay important, but the disk emis- 
sion is still full or partly obscured (very close to the torus 
one can expect absorbing material able to absorb disk emis- 
sion and partly non-disk emission) ; and 3) for a low-inclined 
disk the absorption coming from the torus is negligible, then 
disk and non-disk emission can be fully detected. 

In the case 2), the fraction of the disk emission may be 
too weak. Consequently, it is very hard to extract possible 




Figure 12. Obscuration of the disk emission: l)torus, 2)absorb- 
ing material around the torus and 3) the region without absorp- 
tion (see text in §5). 

disk parameters (see Fig 8, the dots denoted as triangles 
mainly have higher inclinations than 25 degrees). 

5 CONCLUSIONS 

Here we investigated the hidden disk emission in single 
peaked line profiles of type 1 AGN in order to find any 
indication that the disk emission is present in non double 
peaked broad lines. 

Finally, we outline the following c onclusions: 

1) As it was mentioned earlier l|Popovic et all l2004l : 
iBon et al]|2006l ). a two component model (disk -I- non-disk 
region) can well describe the majority of the observed single 
peaked line profiles. 

2) After comparing simulated and observed line profiles, 
there is an indication that the disk emission may be present 
also in single peaked broad line profiles, but it is mainly 
smaller than 50% of the total line flux. 

3) The estimated inclination of the disk indicates a low 
inclined disk, with inclinations i < 25° that might be caused 
by the torus and/or absorbing material around the torus. 

4) In some of the observed Ha line profiles there is a 
blue asymmetry that can be fitted considering a blue-shift 
of the non-disk region of ~800 km/s. This may indicate 
an outfiow in the non-disk region in some AGN from the 
sample. 
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APPENDIX A: THE CASES OF THE 
DIFFERENT GAUSSIAN WIDTHS 

As we mentioned in §2.2, we fixed the widths of Gaussian 
(that represents the n on-disk component) t o a = 1000 km/s. 
Recall the results of I Popovic et al.l (|2004f ) and iBon et al] 
(|2006| ). where the o of the Doppler broadening from the non 
disk component (the Gaussian representing ILR) were esti- 
mated to be from 300 to 1700 km/ s . This is in a good agree- 
ment with the results of IHu et al.l (|2008h . They used a two 
component model assuming that a broad line can be repre- 
sented by two Gaussians (one very broad, and another with 
intermediate velocities, corresponding to ILR and VBLR re- 
spectively) that FWHM's from VBLR and ILR components 
are in correlation, as FWHM(ILR) ~ 0.4 FWHM(VLBR), 
and FWHM (ILR ) ranged from 1 000 to 4000 km/s {a ~ 400 
to ~ 1700 km/s) (iHuet al.ll2008l ). 

But, to have an impression how a of the non-disk com- 
ponent can aflect composite profiles, we performed simula- 
tions using different a. We started from o ~ 400 to ~ 1700 
km/s which corresponds FWHM ~ 1000 to ~ 4000 km/s, 
with a step of A FWHM = 1000 km/s. The disk parame- 
ters are taken as Ri„„=600 Rg, Rotit=4000 Rg, emissivity 
p = 3, and we considered different inclinations and Q. In 
Figs. IA1IIA3I we presented calculated parameters fc2o, ^30 
vs. fcio for Q=0.3, 1 and 2, taking o — 500 km/s, 1000 km/s 
and 1700 km/s and different inclinations. Also, the elliptical 
surface shown in Figs . I A H A3I presents the surface where the 
measured data from the sample are located. As one can see, 
different values of a can affect obtained results for Q and in- 
clination, but it is interesting that measured data indicates 
a ~ 1000 km/s well fit the surface of measured values in 
all three considered cases. Also, for the very broad non-disk 
component (a > 1700 km/s) it is hard to use this model to 
estimate the disk parameters (especially for small inclina- 
tions) . 
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Table 1. Measured k and estimated inclination and ratio of contributions of the non-disk and disk emission. 
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1, 


.92 


1 


.49 


19.05 


0.1 


1.1 


0.1 


SDSS,J1223+0240 


0.0722 


2, 


,17 


1, 


.69 


1 


.41 


14.75 


0.1 


1.3 


0.2 


SDSSJ0755+3911 


0.0335 


2, 


,41 


1, 


.77 


1 


.45 


17.5 


0.05 


1.3 


0.2 


SDSSJ1141-I-0241 


0.0459 


2, 


,78 


1, 


.96 


1 


.52 


19 


0.05 


1.2 


0.01 


SDSSJ1122-I-0117 


0.0394 


2, 


,70 


1, 


.85 


1 


.45 


19 


0.05 


1.25 


0.05 


SDSSJ1243-I-0252 


0.0767 


2. 


,50 


1, 


.81 


1 


.50 


18.85 


0.5 


1.2 


0.2 


SDSSJ0832-I-4614 


0.0605 


3, 


,00 


1, 


.90 


1 


.41 


19 


0.05 


1.25 


0.05 


SDSSJ0840+0333 


0.0525 


2, 


,04 


1, 


.53 


1 


.33 










SDSSJ1510+0058 


0.0359 


2, 


,35 


1, 


.79 


1 


.50 


17.05 


0.3 


1.1 


0.2 



12 E. Bon et al. 



Table 1. continued 



SDSSJOllO-1008 


0.078 


2.39 


1 


82 


1 


49 


17.45 


0.1 


1.15 


0.15 


SDSSJ0142-1008 


0.0303 


2.33 


1 


63 


1 


32 


19 


0.05 


2.4 


0.01 


SDSSJ1519+5208 


0.0693 


2.53 


1 


74 


1 


42 


18.85 


0.1 


1.6 


0.1 


SDSSJ0013-0951 


0.0738 


2.53 


1 


80 


1 


53 


19.1 


0.15 


1.2 


0.3 


SDSSJ1535+5754 


0.0615 


2.12 


1 


64 


1 


36 


15.5 


1.15 


2.5 


0.5 


SDSSJ1654+3925 


0.0419 


2.23 


1 


67 


1 


37 


17.05 


0.9 


2.3 


0.3 


SDSSJ0042-1049 


0.0581 


2.52 


1 


78 


1 


43 


19.15 


0.3 


1.45 


0.15 


SDSSJ2058-0650 


0.0904 


2.50 


1 


81 


1 


38 


19.25 


0.2 


1.65 


0.15 


SDSSJ1300+6139 


0.0522 


2.10 


1 


73 


1 


42 


14.55 


0.2 


0.85 


0.05 


SDSSJ0752+2617 


0.0948 


2.74 


1 


87 


1 


48 


19 


0.05 


1.25 


0.05 


SDSSJ1157+0412 


0.0822 


2.84 


1 


89 


1 


42 


19 


0.05 


1.25 


0.05 


SDSSJ1139+5911 


0.0854 


2.30 


1 


70 


1 


39 


18.55 


1.7 


1.75 


0.15 


SDSSJ1345-0259 


0.0279 


2.00 


1 


53 


1 


30 


20.45 


0.01 


3.7 


0.1 


SDSSJ1118+5803 


0.0613 


2.32 


1 


73 


1 


46 


16.65 


0.1 


1.3 


0.3 


SDSSJ1105+0745 


0.0734 


2.45 


1 


68 


1 


43 


18.45 


0.6 


1.65 


0.35 


SDSSJ1623+4804 


0.0449 


2.33 


1 


71 


1 


37 


20 


0.15 


2 


0.2 


SDSSJ0830+3405 


0.0696 


2.00 


1 


60 


1 


35 


12.1 


0.15 


2.05 


0.15 


SDSSJ1619+4058 


0.0335 


2.12 


1 


62 


1 


34 


17.75 


1.6 


3.1 


0.3 


SDSSJ0857+0528 


0.0379 


2.39 


1 


77 


1 


42 


17.2 


0.05 


1.4 


0.1 


SDSSJ1613+3717 


0.0586 


2.49 


1 


77 


1 


47 


19.05 


0.6 


1.35 


0.25 


SDSSJ1025+5140 


0.0623 


3.00 


2 


11 


1 


54 


19 


0.05 


1.15 


0.05 


SDSSJ1016+4210 


0.0553 


2.63 


1 


83 


1 


42 


19.05 


0.01 


1.35 


0.05 


SDSSJ1128+1023 


0.0504 


2.88 


2 


00 


1 


53 


18.95 


0.01 


1.15 


0.05 


SDSSJ1300+5641 


0.0718 


2.75 


2 


00 


1 


58 


19 


0.05 


1.1 


0.01 


SDSSJ1538+4440 


0.0406 


2.29 


1 


79 


1 


45 


18.2 


2.05 


1.1 


0.1 


SDSSJ1342+5642 


0.0728 


2.46 


1 


88 


1 


54 


18.45 


0.4 


0.95 


0.15 


SDSSJ1344+4416 


0.0547 


2.68 


1 


89 


1 


42 


19 


0.05 


1.25 


0.05 


SDSSJ1554+3238 


0.0483 


2.03 


1 


48 


1 


26 


20.4 


0.05 


3.4 


0.01 




Figure Al. Simulated ^20 vs. feio (left) and fcso vs. fcio (right) for different cr for the non-disk component. The elliptical surface denotes 
the surface of measured values from the sample. The disk parameters are taken as: Rinn=600 Rg, Rout=4000 Rg, (7^= 1000 km/s 



APPENDIX B: THE TWO-COMPONENT 
MODEL VS. TWO-GAUSSIAN FITTING 

As we mentoned above, there are several possibilities for 
the BLR geometry, i.e. for geometry of the VBLR and ILR. 
In principle, there are many analyzes based on two broad 
Gaussian in order to explain phys ics and geometry of the 
BLR (see e.g. iMarziani eral]|2009l) . 

Here we briefly discuss the fitting using the two- 



component model (as it was described in IPopovic et al.l 
120041^ and two Gaussians. As an example here we show 
the best fit with the two-component model (Fig. IBlb .. 
left) and with the two-Gaussian one (Fig. IBlb ). As one 
can see in Fig. IBll both models can well fit the com- 
plex line profiles. The obtained kinematical parameters 
are: 1) two-Gaussian FWHMvsla =3760 km/s (redshifted 
~560 km/s), FWHM/Lfl =1390 km/s (in the center) and 
Q,=P ilr/^vblr 2) two-component model with the 
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disk emission FWHMo^fe =1620 km/s (redshifted ~0 
km/s), FWHM/iH =1360 km/s and Q^F/Lfl/Fjj.^fc «0.8. 
The obtained kinematical parameters in both cases indicate 
that VBLR has higher velocities, but in the case of the as- 
sumed disk emission the random velocity is comparable with 
one present in the ILR. 

Comparing the VBLR Gaussian and disk component we 
found that mainly the VBLR Gaussian should be shifted to 
the red in order to fit complex BEL profiles (similar as i n 
the case of so called Pop B objects, see ISulentic et aLll2009^ ■ 
while the disk component is more consistent with the ILR 
component. 

Additionally, we did several tests taking into account 
the two-Gaussian model, fixing one that represents the ILR 
emission at FWHM=1000 km/s, and changing the width of 
the one representing the VBLR from 1000 km/s to 10000 
km/s taking different Q=Filr/Fvblr~0.3, 0.5, 1, and 2 



(see Fig. IB2p . We also measured ki for such modeled profiles 
and compared them with measured from observed profiles 
(Fig. IB2|I . In Fig. IB2l the lines represent modeled values, and 
crosses measured ones. As it can be seen such model may 
describe majority of the observed profiles, but there is a big 
difference between widths for the VBLR component esti- 
mated using ^20, 30 vs. fcio. Also, in this case we obtain that 
larger fraction of observed AGN has Q=Filr/Fvblr <1, 
i.e. that the VBLR component is dominant in line profiles. 
Comparing these two models (see Figs. 6 and lB2|) . we found 
that the two-component model with the disk emission gives 
more consistent results (Q and i, see Table I) 

In principle for some kind of investigation it is useful 
(and quit simpler) to use the two-Gaussian fit, but here we 
aim to investigate possible presence of the disk emission 
in single-peaked lines, and physically the two-component 
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Figure Bl. The best fit with two broad Gaussian (left) and with two-component model (right) of 3C273 ('see lPopovic et aLll2004f) where 
one component is assumed to be emitted from the disk (a double peaked component below). 
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1.5 2 2.5 3 3.5 ' 1.5 2 2.5 3 3.5 

Figure B2. The measured width ratios (crosses) and simulated values (lines) from the two-Gaussian model for the different contribution 
of the disk emission to the total line flux (Q = Fi lh/ Fy blr=0.3, 0.5, 1, and 2). The width of the ILR component is fixed as 1000 
km/s, and with the VBLR component is taken as 2000 km/s (full circles), 3000 km/s (open circles), 5000 km/s (open triangles) and 
10000 km/s (full triangles). 
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model with a disk emission seems to better explain the na- 
ture of AGN. 
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